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In this paper we discuss decay of superfluid currents in boson lattice systems due 
to quantum tunneling and thermal activation mechanisms. We derive asymptotic 
expressions for the decay rate near the critical current in two regimes, deep in the 
superfluid phase and close to the superfluid-Mott insulator transition. The broad- 
ening of the transition at the critical current due to these decay mechanisms is more 
pronounced at lower dimensions. We also find that the crossover temperature below 
which quantum decay dominates is experimentally accessible in most cases. Finally, 
we discuss the dynamics of the current decay and point out the difference between 
low and high currents. 

I. INTRODUCTION 

Some of the most intriguing questions in low temperature physics concern the ways 
in which superconductors lose their superconducting properties, because of thermal or 
quantum fluctuations. Mike Tinkham has long been fascinated with these issues, and has 
done much to advance our understanding of the subject. 

An early contribution in this area was the work of Newbower, Beasley and Tinkham 



on effects of fluctuations on the superconducting transition of tin whisker crystals^. Ex- 
perimental data were compared with theories of thermally activated phase slips, both 
in the linear regime and in the non-linear regime of finite current flows. More recently, 
Tinkham and collaborators studied the loss of superconductivity in very thin wires of 
MoGe, deposited on carbon nanotubes, where quantum fluctuations are involved^. Re- 
lated work from Tinkham's laboratory, in recent years, has elucidated the breakdown of 
superconductivity in ultrasmall metallic grains, measured by the even-odd alternation of 
Coulomb-blockade energies^ , vortex motion and resistance in high temperature super- 
conductors^ , and critical currents in frustrated arrays of Josephson junctions.—. 

The decay of supercurrents in liquid 4 He and in Bose-Einstein condensates of ultracold 
atoms has much in common with the decay of superconductivity. Concepts of flux-line 
motion, and of phase slips due to thermal or quantum fluctuations, appear in both cases. A 
new dimension has been added to the subject by recent experimental advances, where cold 
atoms have been trapped in a region that contains a spatially periodic potential, produced 



by optical standing waves (see for example 



lCj). The ability to vary continuously the 



parameters of the system, by changing the strength of the periodic potential, as well as 
by varying the number of trapped atoms and the shape of the overall confining potential, 
allows one to explore new regimes of parameters and to make more precise confrontations 
between theory and experiment. In turn, these developments give added urgency to the 
theoretical study of supercurrent decay. 

In the present paper, we discuss similarities and differences between the decay of su- 
percurrents in superconductors and systems of trapped atoms, and we present some new 
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results for the latter. Specifically we consider certain experimental procedures which have 
become standard in systems of ultra cold atoms. In the first scheme, a condensate is pre- 
pared on a lattice with a specified intensity, when the lattice is suddenly accelerated to a 
finite velocity. In other words, a moving condensate is prepared in the lattice frame, essen- 
tially fixing the phase gradient. A similar experiment in superconductors would involve 
threading a flux through a closed superconducting loop. Such sudden lattice boosts were 
applied by several groups to demonstrate a dynamical instability of the superfluid when 
the imposed phase gradient exceeds ir/2 per unit celt^. A related experiment involves 
tilting the lattice, thus subjecting the atoms to a linear potential. This is equivalent to 
imposing a constant voltage on a superconductor. The technique was used, for example, 
to demonstrate Bloch oscillations of a condensate 11 . In a third experimental sequence, 
one can prepare a moving condensate, then continuously increase the depth of the lattice 
toward the superfluid- Mott insulator transition. 

The response of the atomic system to the perturbations, can be measured, by direct 
observation of the time evolution. Decay of the current, for example, is observed by re- 
peated experiments, where atoms are released from the trap after varying waiting periods. 
The phase gradient in the superfluid at the time of release may be inferred from a time 
of flight measurement of the momentum distribution. This should be contrasted with 
superconductors, where measurements probe I-V characteristics. 

Besides the differences in the experimental observation procedures, there are unique 
features of trapped atom systems which influence the physics of super current decay. First, 
to a very good approximation such systems can be considered perfectly clean. Super- 
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currents decay only due to breaking of Galilean invariance by the periodic potential. 

A second feature that distinguishes the dynamics of ultra cold atoms is their nearly 
perfect isolation from the environment. Strictly speaking they are always under-damped. 
However, we are usually interested in the dynamics of a subset of the system degrees 
of freedom, such as the super-current. How much the dynamics of the the interesting 
variables is damped, depends solely on the remaining system degrees of freedom rather 
than on external dissipation sources. In superconductors on the other hand, the order 
parameter is invariably in strong contact with other degrees of freedom, and therefore its 
dynamics is almost always overdamped. 

II. CRITICAL CURRENT IN THE SUPERFLUID PHASE 

Ultra cold atoms in an optical lattice, confined to the lowest Bloch band are described 
by the well known Bose Hubbard Hamiltonan: 

H — — J J2( a l a i + H.c.) + ni{ni ~ X )' W 

(ij) i 

where J and U are the hopping amplitude and the on-cite repulsive interaction, (ij) 
denotes pairs of nearest neighbor sites. Another implicit parameter in this Hamiltonian is 
the average number of bosons per site, N. In this paper we shall be primarily concerned 
with the case where A is a large integer. We shall address to separate regimes, the first is 
defined by the conditions UN 2 ^> JN ^> U, while the second regime corresponds to the 
superfluid near the transition to a Mott insulator (UN 2 ^> JN ~ U). 

If the condition UN ^> J holds, then the interactions are sufficiently strong to suppress 
amplitude fluctuations of the order parameter, and (JIJ can be mapped to the quantum 
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phase model: 




(2) 



The additional condition JN ^> U ensures that the system is far from the superfluid- 
insulator transition, and facilitates a semiclassical approximation. In the classical limit 
the boson creation and annihilation operators can be treated as complex numbers subject 
to discrete Gross-Pitaevskii equations^: 



where the set O contains the nearest neighbors of site j. In the quantum rotor limit 
UN ^ J the number fluctuations can be integrated out leaving us with only the equations 
of motion for the phase <ftj = arg-i/^: 



Alternatively Eq. immediately follows from the Hamiltonian (J2J). Both equations 
(jUJ) and (HJ) can support stationary current carrying states ipj oc exp(ipxj). A simple 
linear stability analysis of show n 14 ' 15 that these states become unstable towards small 
perturbations when the phase twist exceeds tt/2 per unit cell. The onset of this instability 
is signaled by appearance of imaginary frequencies. So 7r/2 is the critical phase twist above 
which a uniform superfluid state breaks down. This instability was recently observed 
experimentally^. 

In principle, one can identify another type of instability, characterized by appearance 
of negative frequencies, in systems described by Eq. (jHJ^ 4 -. In general this occurs at a 
phase twist p* < tt/2. However, in the quantum rotor limit UN ^> J, where we work, the 
two instabilities coincide. 




(3) 




(4) 
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While the modulational instability occurs precisely at p = tt/2 for JN » U, we 
expect that the current can decay at smaller momenta due to either quantum or thermal 



fluctuations (see also Ref. Ilfil). We envision the following experimental scheme to observe 
this. The condensate is either boosted to a state with a certain phase gradient or gradually 
accelerated. Following the boost or while the system is accelerating we probe the evolution 
of the phase gradient. If the system is sufficiently close to the modulational instability, 
i.e. p is slightly below tt/2, the coherent motion of the condensate is expected to decay, 
the larger the phase gradient, the faster this decay will occur. 

The other regime we shall address, is that of the superfluid close to the quantum phase 
transition to a Mott insulator at commensurate filling (i.e. JN ~ U). For simplicity 
we also assume that N ^> 1. Since the coherence length £ diverges at the transition, 
one can use a continuum semiclassical description, after coarse-graining the system, to 
describe the static and dynamic properties of the condensate. At commensurate filling the 
appropriate quantum action written in terms of the superfluid order parameter reads : 17 i 18 



S = C d d xdt 



dtp 2 



IWf + H 2 -^I 4 , (5) 



dt 

where length is measured in units of £ and time in units of £/c, with c the sound velocity. 
C is a numerical prefactor. The bare parameters £, c, and C can be found using a 
mean-field approximation^. For the cubic d-dimensional lattice they read: 

j= t 1 =, c = 2 JNV2d, C = . ^ (1 -u)^ , (6) 
y/2d(l-u) 2(2^/2 ^ ' ' y ' 

where we introduced the dimensionless interaction u = U/U c with U c = 8JNd the critical 
interaction strength in the mean field approximation. The action (j3J) correctly describes 
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low energy dynamics of the system in the vicinity of the phase transition, only if the 
couplings £, c, U c and C are properly renormalized. While in three dimensions the effects 
of such renormalization should be weak, in two and especially one dimensional cases they 
strongly modify the couplings and the critical exponents. The bare mean field parameters 
can then be used only as an estimate. Note that the dimensionless part of the action (jSJ) 
is general, and so are the conclusions we reach in this paper, once the renormalized, 
rather than mean field parameters are used in the prefactors. The action (jHJ) is obviously 
extremized by stationary current-carrying states: ifi p (x) = yl — (j»£) 2 e tpx ^- It is easy to 
checks that these states are stable with respect to small fluctuations for p < p c = l/(£\/3). 
Since £ diverges at the phase transition, the critical phase twist vanishes at that point as 
it should. 

A possible experimental procedure to measure the decay rate at low currents follows. A 
condensate with a specified phase gradient is prepared in a weak lattice (small u). Then, 
the lattice potential is gradually increased in time, driving the system closer to the Mott 
phase. This, in turn, results in the increase of the correlation length £ and in decrease of 
the critical momentum p c . As p c approaches p the superfluid current is expected to decay 
either due to quantum or thermal fluctuations. 

III. DECAY OF THE SUPERFLUID CURRENT 

In this section we describe how the superfluid current decays in a lattice when p is 
below p c . We shall consider first the Gross-Pitaevskii regime JN >> U in the quantum 
rotor limit UN » J . due to either thermal or quantum effects. 
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A. Gross-Pitaevskii regime 

1. Quantum Decay 

We first discuss the decay of superfluid currents within the quantum phase model, well 
away from the superfluid insulator transition (JN ^> U). The action corresponding to 
the quantum rotor model O) is given by 



5 = ' dT E (v 1 ) 2 - E 2^008^ - <P f ), (7) 



2U V dr 



or after the rescaling r — > rj yUNJ : 



v JN 

where 

: cos(0j — (9) 



j v 7 0,i'> 



To leading order in v/C// JN, which plays the role of the effective Planck's constant for 
this problem 1 ^ the tunneling rate depends on the action Sb, associated with the bounce 
solution of the classical equations of motion in the inverted potential 2 ^: 

T oc e~ s \ (10) 

Clearly the action should vanish at p = tt/2, since at this point the spectrum becomes 
unstable and the tunneling barrier disappears. Deep in the superfluid regime U / JN <^ 1, 
the tunneling is effective only if p is close to tt/2, where the product W JN/Us is not too 
large. In this case one can make further progress in calculating the tunneling action by 



expanding (JUJ) up to cubic terms in phase differences 



3 • 



s « 2J / ^ r 



1 f dcf) \ ^ 1 

(11) 



Here we explicitly split the site index into longitudinal (J) and transverse (k) components. 
Also, for the convenience, we shifted the phase c/^k — ► (ftj^+pj so that the metsatbale state 
now corresponds to (pj^ = 0. Note that at p — > 7r/2 only longitudinal modes become soft, 
acquiring a prefactor cosp in front of the quadratic term in the action. This implies that 
we can safely use continuum approximation for the phases along the transverse directions. 
Then instead of ([lip we derive: 



(12) 



In this equation x denotes transverse coordinates which reside in a d — 1 dimensional 
space. Upon rescaling 

= cos(p)4>, T = J-= , x = ^z^L, (13) 
y/cos(jp) V cos P 

the action f!12l) simplifies further: 

s ^ (n/2-p)^§ d , (14) 

where 



d-l 

§d = 2— 



2 



(15) 



is just a number, which is determined only by the dimensionality d. We will pro- 
vide its detailed variational derivation elsewhere 21 and here just quote the results: 
s\ ~ 7, §2 ~ 25, §3 ~ 90. ^From the scaling (fT3*j) it is obvious that the characteris- 
tic transverse dimension of the instanton x scales as (it/2 — p)~ x l 2 1, justifying the 
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continuum approximation. Above d = 6 the tunneling action would experience a dis- 
continuous jump at p = tt/2. However, since we deal with d < 3, the action always 
continuously vanishes at p — > tt/2. In this way we derive the asymptotic decay rate of a 
uniform current state near the modulation instability: 

T oc e" 5d (16) 

with 

S d ^~s d ^^/2-p)^. (17) 
2. Thermal Decay 

To calculate the exponent, characterizing the thermal decay rate one has to compute 
the difference of energies of the metastable state and the saddle-point, which separates 
two adjacent metastable minim a 22 i 23 i 24 . Both saddle-point and metastable configurations 
are the stationary solutions of the equations of motion (J3J): 

^sin^-^) = 0. (18) 
fceo 

The metastable state corresponds to the uniform phase twist: 4>j,\a = 3P- The saddle-point 
state relevant for the durrent decay can be easily found in one dimension: 

{p'j j < 

(19) 
TT+p'U-2) j>\ 

where p' ~ p — (n — 2p)/M if we use periodic boundary conditions for the system with M 
sites. The energy difference between the two states in the limit M — > oo is 

AE = 2JN(2cosp - sinp(7r - 2p)). (20) 
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Correspondingly the decay rate is proportional to 

Y e -/3A£ _ e -2JAr/3(2cosp-(7r-2p)sinp)_ ^1) 

In particular when p — > n/2 we have: 

r(xe -|iV W 2-p)3_ (22) 

In higher dimensions we cannot find the energy of the saddle-point explicitly for all values 
of p. However, similarly to the quantum case, at p close to 7r/2 we can expand the energy 
functional up to cubic terms in phase differences 



(23) 



where <frj(x) is the nontrivial solution of the corresponding Euler-Lagrange equations van- 
ishing at x — > oo. We again shifted the phase: 4>j ~~ * 4>j +P3-, so that <f>j = corresponds 
to the current state. After rescaling tpj = cos(p) <j)j and x = x\f2j ^/cosp we find: 



E d ^ JN2^(p c -p) 1 ^J2 J dd ~ l3c 



1,7 

2 \ dx I ' 1 3 



.(24) 



Note that the integral in the expression above coincides with Sd-i up to a number 2 2 . 
So we immediately conclude that 

E d « s d _i JNy/2(n/2 - p) ^ . (25) 

Note that the activation energy characterizing the thermal decay vanishes faster as p — > 
7r/2 than the tunneling action. It implies that thermal fluctuations become increasingly 
important and dominate the decay of superfluid current as the system approaches the 
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dynamical instability. Comparing the ratio E^/T and the tunneling action (|17j) we obtain 
the crossover temperature: 

T* « £±l Cv /7r/2-p (26) 

at which the quantum and thermal decay rates coincide. Here c = >/2{7 JN is the sound 
speed in equilibrium (i.e. p — 0). Alternatively, we can fix the temperature to obtain the 
momentum crossover scale p* at which thermal and quantum decay rates coincide: 



i/2 - p *n^j {-) ■ (27) 

At phase gradients larger than p*, thermal decay dominates. The tunneling action (fTTj) 
at this value of momentum is given by 



Jd-iJ V U \c 

If S*i 3> 1, then at the crossover momentum the current decay is exponentially suppressed 
and will be nonzero only at p closer to ir/2. Then the thermally activated phase slips will 
dominate the decay process and quantum tunneling can be ignored. In the opposite limit 
5,5 <C 1 the current will decay at p < p* due to quantum process and the temperature 
effects are unimportant. The characteristic crossover temperature separating quantum 
and thermal decay regimes is thus: 

/ U \ 2(6^d) 

T q - Ac~s d - (-) , (29) 

where A is a numerical constant of the order of one. Note that for all relevant dimensions 
d < 3 the last multiplier is always of the order of one because of the small exponent 
1/(12 — 2d). Therefore, the crossover temperature T q is of the order of the sound velocity 
(or equivalently Josephson energy). 
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IV. CURRENT DECAY IN THE VICINITY OF THE MOTT TRANSITION 

Let us now address decay of supercurrents in the regime where JN ~ U and large 
integer filling N. As we already argued, in the vicinity of the Mott-insulator phase 
transition the correlation length £ becomes large compared to the lattice constant. One 
can therefore use a continuum description of the problem (J3J). The relativistic dynamics of 
(0) is a special feature of the commensurate transition. Ordinary superfluids are described 
by a similar action, but with a kinetic term including first time derivatives. 

The Euler Lagrange-equations derived from © admit stationary current carrying so- 
lutions of the form 



We emphasize again that x is measured in units of the correlation length £. The current 
state becomes unstable at p > p c = l/>/3£. Below p c the current can still decay due to 
quantum or thermal fluctuations. 

For the thermal decay only the static part of the action needs to be considered. Then 
there is no difference between our problem and the current decay in ordinary superfluids 
described by the Ginzburg-Landau free energy. In particular, in context of supercon- 
ducting wires, the exponent characterizing the current decay rate in one dimension was 
computed by Langer and Ambegaokar— , and the prefactor setting the time scale was later 
found by McCumber and Halperin 24 . In three dimensions the asymptotic behavior of the 
corresponding exponent at p —>■ was obtained by Langer and Fished. However, here 
we are interested in the opposite limit p — ► p c . 

For both the thermal and the quantum cases we will use the scaling approach sue- 




(30) 
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cessfully applied above for the quantum phase model in the Gross-Pitaevskii regime. We 
expand the action to cubic order in the amplitude (rj) and phase (0) fluctuations, about the 
metastable minimum, and integrate out the gapped amplitude mode. After the rescaling: 



x z 3 3 / 4 



x 



VtVK^p (3i) 



2 3 1 /V£ v ^=p' Q^Pc-pV 
the action to leading order in p c — p becomes: 

o9/4-d r 

S = C^^e /2 ' d (Pc-p) 5/2 ~ d J dzdx(V0) 2 +(^0) 2 -(40) 3 « A d {l-u) l ! A (p c -pf' 2 - d 

(32) 

where z denotes all the transverse coordinates relative to the current direction, including 
time. V = (<9 Z , d x ) is the gradient in d + 1 dimensions. Accordingly, the quantum decay 
rate is given by Tq oc exp(— Ad(l — u) 1 / A {p c — p) 5 / 2 ~ d ). A variational calculation 2 ^, yields 
A\ pa 18.4 and A2 pa 8.4. As before, to calculate the thermal decay rate one simply has 
to substitute d — > d — 1, so that 

T T (d) oc exp (- J -^L(2d)-*l\l - ufl'A^ipc - py' 2 - d \ . (33) 



In the one-dimensional case the relevant constant Aq pa 12.56. It is interesting to contrast 
these results with the asymptotic decay rate ()17|) . found in the Gross-Pitaevskii regime. 
First we observe that the tunnelling action close to the Mott insulator vanishes as a 
smaller power of p c — p. Moreover, for d = 3 the scaling hypothesis for the quantum 
decay rate breaks down, suggesting that S is discontinuous at the critical current and 
is dominated by the fluctuations of a finite (not diverging as p — > p c ) length scale. We 
therefore expect that in three dimensions at zero temperature the instability marks a 
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sharp localization transition. At finite T, thermal fluctuations broaden this transition, 
because the activation barrier energy vanishes at p c for d < 7/2. 

The quantum to thermal crossover for a given dimensionless interaction and phase 
gradient is found by comparing the two decay rates. In one and two dimensions we find 

r(p) =(2^^- p) ' (34) 

In three dimensions T* = because the quantum decay is effectively suppressed. As 
discussed above for the quantum phase model, there is a more useful p independent, 
crossover temperature scale. Using the same arguments as in the Gross-Pitaevskii we can 
find the temperature separating the quantum and thermal decay regimes: 

7-2d 1 

T q ~ JNA d ^A d 5 - 2d (l - u)-™=a. (35) 

We see that near the Mott transition the crossover temperature strongly depends on 
interaction u. Thus as u — > 1 then in one and two dimensions T q — > oo so the quantum 
decay always dominates over the thermal. In particular in two dimensions we find T q ~ 
0.03 JN / y/1 — u and in the one dimensional case T q ~ 0.1JA^/(1-m) 1 / 6 , i.e. the crossover 
temperature is very high and the thermal decay is unimportant. 

V. DYNAMICS OF THE CURRENT DECAY 

Except for the zero temperature decay in the continuum relativistic model in three 
dimensions, there is no sharp transition between superfluid current carrying state and the 
insulating state with no current. Indeed, in all other cases the tunneling action and the 
energy barrier continously vanish as the system approaches the point of the modulation 
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instability. Thus, instead of a sharp transition boundary we can define a broad crossover 
region, say 1 < Sd < 3, which separates the superfluid phase with a relatively slow current 
decay and the insulating phase with a fast decay. The fact that the transition is broad 
does not imply, however, that within a single experiment a gradual current decay will 
be detected as the system slowly tuned through the crossover region. The tunneling and 
thermal decay rate define a probability of creating a single phase slip per lattice site. In 
the following evolution the phase slip can either rapidly release its energy into phonon 
(Bogoliubov's) modes and bring the system to a next metastable minimum with a lower 
current, or it can trigger the current decay in the whole system. 

In a closed system, i.e. with no coupling to the environment, these two regimes are well 
defined because the damping of the phase slip comes from the internal degrees of freedom, 
which are completely described by the equations of motion. Also near the critical current 
there should be very little difference between the overdamped and underdamped regimes 
in the thermal and quantum cases. This is because the energy barrier is very small and 
the classically allowed motion after the tunneling event starts very close to the metsatable 
maximum. To see which of the regimes is realized in our systems we numerically solve 
the classical (Gross-Pitaveskii) equations of motion We start from a uniform current 
state in a periodic lattice. To allow for a current decay we add small fluctuations to 
the initial values of the classical fields ipj(t = 0). This mimics the effect of thermal 
fluctuations. In Fig. [T] we plot the computed current versus time for a one-dimensional 
array of M — 200 sites. Initially the system is asumed to be noninteracting (U = 0) 
macroscopically occupied state with a given phase gradient p (specifically we consider 
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p = 27r/5 and p = tt/10) and unit hopping. Then interaction is gradually increasing in 
time and reaches some constant value and the we follow the time evolution of the current. 
It is clear from the figure that the phase slips in the smaller current case (p = 7r/10) are 
overdamped leading to gradual decay. With the used parameters the total winding number 
is 2007r/10/27r = 10 and the figure clearly shows that each phase slip decreases the current 
by roughly ten percent. On the other hand for the larger current (p = 27r/5) a single phase 
slip generates immediate current decay in the whole sample and this corresponds to the 
underdamped regime. We will not attempt here to find the precise boundary between 
the scenarios, since it is not the purpose of this paper. We would like to stress that 
near the instability the system is always in the underdamped regime. We checked that a 
similar overdamped to underdamped crossover occurs in other spatial dimensions. So if 
p is not too small, the in a given experiment one will always see a sharp transition from 
the superfluid to the insulating regime. However the precise point, where the current 
decays will depend on the details of the experiment, for example on the rate of change of 
external parameters like tunneling, interaction, or the phase gradient p if the system is 
accelerated. On the other hand in the absence of any fluctuations the transition is very 
sharp and always occurs at p = ir/2. 

We can do a similar numerical analysis at small currents close to the Mott transition 
using the action We find that even close to the critical current the decay is over- 
damped, i.e. the current decays in steps. This result is natural because the size of the 
phase slip in this case is large and hence it is very easy to dissipate the released energy 
into phonon modes. 
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We would like to mention that if the motion of phase slips is under damped then in 
a truly infinite system the current state is always unstable. Indeed the probability of 
a phase slip in the whole system is proportional to its size M. If this single phase slip 
causes triggers the current decay in the whole sample, then obviously a state with a 
uniform phase gradient cannot exist. However, in finite size systems these effects are 
not so crucial, because the decay probability depends exponentially on the couplings and 
current but only linearly in a system size. 

VI. CONCLUSIONS 

The modulational instability can be observed either by increasing momentum of the 
superfluid accelerating the condensate or by increasing lattice potential and driving the 
system closer to the Mott transition, while the condensate is in motion. We showed 
that because of quantum or thermal effects the current decays before the system becomes 
classically unstable. Therefore instead of a sharp transition, there is a crossover region 
where the decay rate grows from being exponentially small to large. The crossover region 
becomes narrower either deep in the superfluid regime (i.e. JN ^> U) or in higher 
dimensions. In particular, in the three-dimensional case we always expect a very sharp 
boundary separating the regions with very weak and very strong decay rates. 

We found that deep in the superfluid regime the crossover temperature separating 
quantum and thermal decay is of the order of the Josephson plasma frequency in all 
dimensions. At small currents, close to the Mott phase, the decay occurs through ther- 
mal fluctuations in three dimensions and through quantum tunneling in one and two 
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dimensions. 

We argue that both overdamped and underdamped dynamics of the current decay can 
be realized in these systems. The underdamped regime corresponds to high currents close 
to p = it/2, while at low currents the dynamics is overdamped. 
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FIG. 1: Current (scaled to one at t = 0) versus time for a one dimensional periodic array 
of 200 sites with two different initial phase gradianets. The evolution is determined solving 
equations of motion © with constant hopping amplitude J = 1 and interaction increasing in 
time U = 0.01 tanh O.Olt for p = 2ir/5 and U = tanhO.Olt for p = tt/10. To get the current 
decay we add small fluctuations to the initial values of the classical fields ipj(t = 0). 



